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1 Introduction 

This goal of this meeting is the search for exotic states of matter in compact as- 
trophysical objects. Over the years many possible signatures for phase transitions 
in compact stars have been suggested, for example unusual masses, radii or cooling 
histories, or sudden changes in the spin frequency. In this contribution we will take a 
very conservative approach and study the predictions of weak coupling QCD for the 
equilibrium and transport properties of the densest phases of QCD. This includes the 
color- flavor- locked (CFL) phase [Ij, and phases that arise from modifications of the 
basic CFL pairing pattern due to the effects of the non-zero strange quark mass. A 
flavor rotation of the CFL condensate leads to kaon condensation (CFL-K) [2], and a 
spatial modulation of the CFL state leads to the meson supercurrent state (curCFL) 

EH]. 

We refer the reader to our recent review ^ for a detailed discussion of quark 
matter phases at lower density. Some of these phases are stable in weak couphng, 
and their properties are rigorously computable in perturbative QCD. This includes 
crystalline color superconductivity [6l [7] , and single flavor spin-one color supercon- 
ductivity [H [9]. Other phases require strong coupling, for example the 2SC phase 
[ini [11], chiral density waves (and the quarkyonic phase) [121 IISI [H], or gapless color 
superconductivity jT5| [T6]. 

The goal of our research program is to compute the properties of all these phases, 
and to determine whether hybrid neutron/quark matter stars containing quark mat- 
ter in the CFL phase or one of the lower density quark phases are consistent with 
observation. Current data on the masses and radii of compact stars are consistent 
with pure neutron stars as well as with hybrid quark matter stars [TTl [18]. There are 
some results that indicate that other observational properties also do not clearly dis- 
tinguish neutron stars from hybrid quark stars. In that case evidence for the existence 
of a high density phase in compact stars may have to come from careful studies of the 
mass-radius relation of compact stars [T9], or from the observation of gravitational 
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waves from binary compact star mergers. Laboratory experiments involving heavy 
ions have estabhshed meaningful constraints on the equation of state of dense matter 
[20] , and the next generation of these experiments may demonstrate the existence of 
a first order transition at densities that are achieved in compact stars. 



Calculations based on weak-coupling QCD indicate that the ground state of three 
flavor baryonic matter at the highest densities is the color-flavor-locked (CFL) phase 
[U HH |22] . The CFL phase is characterized by a pair condensate 



This condensate leads to a gap in the excitation spectrum of all fermions and com- 
pletely screens the gluonic interaction. Both the chiral SU{3)l x SU{3)ji and color 
SU{3) symmetry are broken, but a vector-like SU{3) flavor symmetry remains un- 
broken. 

The gap in the fermion spectrum can be computed in perturbative QCD. The nine 
quarks species (three flavors and three colors) form an octet and a singlet under the 
unbroken SU{3) flavor symmetry. The octet and singlet quarks have gaps Aq and 
2Ao, respectively, where [231IM l25l 126] 



Here, fi is the chemical potential and g is the strong coupling constant. At this order 
in g we can compute the coupling by evaluating the one-loop running coupling at 
the scale /i. At large /x the coupling constant is small and the gap is exponentially 
small compared to the Fermi energy Ef = fi. For densities relevant to neutron stars 
fi < 500 MeV and the coupling is not small. In this regime higher order correction to 
equ. ([2]) are not small and the gap is quite uncertain. It is remarkable, however, that 
both extrapolations of the weak coupling result to the regime of moderate coupling 
[27] as well as model calculations based on Nambu-Jona Lasinio (NJL) or similar 
interactions [28j give gaps in the range Aq — (50 — 100) MeV at a density p ~ 5po, 
where po is the nuclear matter saturation density. 

Perturbation theory can also be used to compute the gluon screening mass. Screen- 
ing arises from the Meissner effect, as it does in ordinary superconductors. We find 



2 Matter at the highest densities 



2.1 The CFL Phase 
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which shows that the gluon screening length is short compared to the coherence 
length ^ ~ Aq^. We conclude that color superconductivity is type I. We note that 
electromagnetism is unbroken, and the photon is not screened. The physical photon 
in the CFL phase is a mixture of the bare photon and the bare gluon, and interesting 
effects occur at the interface of normal nuclear matter and the CFL phase |31j . 



2.2 Effective theory of the CFL phase 

For excitation energies smaller than the gap the only relevant degrees of freedom are 
the Goldstone modes associated with the breaking of SU (3) l x SU{3)r chiral symme- 
try, the U{1)b phase symmetry associated with baryon number, and the approximate 
U{1)a phase symmetry. The interaction of the low energy modes is described by an 
effective field theory. The structure of the effective lagrangian is determined by the 
symmetries of the CFL phase, and the coefficients that appear in the lagrangian can 
be computed in perturbative QCD. We will see that, with one notable exception, 
these coefficients agree with estimates based on naive dimensional analysis. 

The effective lagrangian for the low energy modes has two important applications. 
First, it determines the spectrum and the interactions of quasi-particles. Based on 
this knowledge we can compute the specific heat and the transport properties of the 
CFL phase. Second, the effective theory determines the response of the CFL phase 
to non-zero quark masses and to external fields, such as lepton chemical potentials or 
magnetic fields. The effective lagrangian therefore determines the phase structure at 
non-asymptotic densities. 

The Goldstone mode associated with superfluidity is related to the phase ip of the 
order parameter 

(CC^75<) = e''^%Be''^<l>^ , (4) 

where 0^ parametrizes the color-flavor orientation of the order parameter. The field 
(fi transforms as <y9 — )■ (y9 + a under U{1)b transformation of the quark fields ip — ?• 
exp{ia)tp. At leading order in the weak coupling limit the effective lagrangian is 
completely fixed by Lorentz invariance and U{1)b symmetry. One can show that 
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{do^-f^y-iv^y +..., (5) 



where . . . denotes terms that are higher order in ^f, or terms of the form d^(p^ with 
n > m. In the low energy limit we can expand the lagrangian in powers of dip. We 
also rescale the field as = {3^/TT)(f in order to make it canonically normalized. We 
will refer to as the phonon field. The phonon lagrangian is 

^ = ^(^00)' - lv'id,<py - ^5o0(9,09^0) + Y^(5.</'5"</>)' + • • • , (6) 
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where v = l/y/S is the speed of sound. We observe that three and four phonon 
vertices are suppressed by powers of |(90|//i^. 

The effective lagrangian for the Goldstone modes associated with chiral symmetry 
breaking has the same structure as the chiral lagrangian at T = = 0. The main 
difference is that Lorentz invariance is broken and only rotational invariance is a good 
symmetry. The effective lagrangian is given by 



^Tr 
4 



VoSVoE^ 



+ 



AiTr(MS^)Tr(MSt) + A2Tr(MS^MSt) 



+A3Tr(MSt)Tr(M^S) + h. 



c. 



+ 



(7) 



Here S = exp{i(f)°'X°' / f.^) is the chiral field, /^r is the pion decay constant and M is 
the mass matrix. The chiral field and the mass matrix transform as S — t- LT^R^ and 
M — )■ LMB) under chiral transformations (L, R) G 5^7(3)^ x SU{?>)r. In order to 
determine the structure of the effective theory we treat M as a field, but in practice 
we are interested in the case M = diag(m„, m^, m^). 

The coefficients f.^, v^,, B, Ai, . . . can be computed in weak coupling perturbation 
theory. In the case of fj^ and this is most easily done my matching the screening 
masses for flavored gauged fields. The coefficients B, Ai are related to the quark mass 
dependence of the condensation energy in the CFL phase. 

At leading order in as the Goldstone boson decay constant and velocity are [29] 
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The coefficient B is related to instanton effects. We find E 
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which shows that B is strongly suppressed at large chemical potential. The Ai terms 
receive perturbative contributions and are given by [35] 
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A. = 0. 
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Finally, one can show that Xl = MMy{2pF) and Xr = M^M/{2pf) act as effec- 
tive chemical potentials for left and right-handed fermion. These effective chemical 
potentials appear in the time derivative of the chiral field [2J , 
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Figure 1: Phase structure of CFL matter as a function of the effective chemical 
potential fis = ^I/C^Pf) and the lepton chemical potential /ig, from Kaplan and 
Reddy |36] . A typical value of fis in a neutron star is 10 MeV. 



3 Matter at non-asymptotic density 

At non-asymptotic densities we can not rely on perturbative QCD calculations to 
determine the magnitude of the gap parameter, but the argument that pairing oc- 
curs and that the CFL phase is energetically favored is quite general, and does not 
depend on details of the interaction. The dominant stress on the CFL phase at non- 
asymptotic densities arises from flavor symmetry breaking due to the quark masses. 
We will focus on the physically relevant case rris ^ rrid — irtu- In this case the main 
expansion parameter is m^/(yuA), which is the ratio of the mass correction to the 
Fermi energy of the strange quark over the magnitude of the gap. 



3.1 Kaon condensation 



Using the chiral effective lagrangian we can determine the dependence of the order 
parameter on the quark masses. The effective potential for the order parameter is 



eff 



J TT 



Tr 



Tr(MS"f)) -Tr((MS"f) 



(12) 



The first term contains the effective chemical potential /i<j = m1/{2pp) and favors 
states with a deficit of strange quarks (with strangeness S = —1). The second 
term favors the neutral ground state S = 1. The lightest excitation with positive 
strangeness is the meson. We therefore consider the ansatz E = exp(iQ;A4) which 
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allows the order parameter to rotate in the direction. The vacuum energy is 

V{a) = -fl (^4^) ' + «)'(cos(«) - 1) j , (13) 

where {'wPkY ~ i^^i/ fn)^{^ + ^s)- Minimizing the vacuum energy we obtain 

{1 /is < m% 

Using the perturbative result for Ai we can get an estimate of the critical strange 
quark mass. We find 

msicrit) = 3.03 ■ m]/^A'^^^, (15) 

from which we obtain ms{crit) ~ 70 MeV for A ~ 50 MeV. This result suggests 
that strange quark matter at densities that can be achieved in neutron stars is kaon 
condensed. The phase structure as a function of the strange quark mass and non- 
zero lepton chemical potentials was studied by Kaplan and Reddy [36], see Fig. [H 
We observe that if the lepton chemical potential is non-zero charged kaon and pion 
condensates are also possible. 



3.2 Fermions in the CFL phase 

Kaon condensation occurs for /i^/A ~ y^mrnj//! ^ 1. For conditions relevant to 
neutron stars /i<j/A can get significantly larger, reaching /i^/A ~ 1. In this case some 
of the fermion modes may become gapless or almost gapless ^16]. In order to study 
this regime we have to include fermions in the effective field theory. The effective 
lagrangian for fermions in the CFL phase is ES] 

£ = Tr (n^v^D^n) - Dli (n^^"^^ {A^, N}) - FTr (n^v^'-i^ [^^, A^]) 

-I- ^ {(Tr {NlNl) - [Tr (iVi)]') - {L ^ R) + h.c] . (16) 

Nl^r are left and right handed baryon fields in the adjoint representation of flavor 
5'[/(3). The baryon fields originate from quark-hadron complementarity [39]. We 
can think of as describing a quark which is surrounded by a diquark cloud, A'^^, ~ 
The covariant derivative of the nucleon field is given by -D^jA^ = d/^^N + 
i[Vfj^,N]. The vector and axial- vector currents are 

V. = ~ {Cd,e + eo,^} , A, = -'-i (V,st) ^, (17) 

where ^ is defined by .^^ = S. F and D are low energy constants that determine the 
baryon axial coupling. In perturbative QCD we find D = F = 1/2. The effective 
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Figure 2: This figure shows the fermion spectrum in the CFL phase. For = there 
are eight fermions with gap A and one fermion with gap 2 A (not shown). Without 
kaon condensation gapless fermion modes appear at /i^ = A (dashed hues). With 
kaon condensation gapless modes appear at fig = 4 A/3. 

chemical potentials {X^jXr) appear in the covariant derivative of the nucleon field. 
We have 



where Xl = MMy{2pF) and Xr = M^M/{2pf) as before. {Xl,Xr) covariant 
derivatives also appears in the axial vector current given in equ. ( !T7|1 . 

We can now study how the fermion spectrum depends on the quark mass. Since the 
field N has the quark numbers of the baryon octet and singlet we will use (p, n, E, S, A) 
to label the fields. In the CFL state we have ^ = 1. For /ig = the octet has an energy 
gap A and the singlet has gap 2A. As a function of fig the excitation energy of the 
proton and neutron is lowered, Up^n = A — /i^, while the energy of the cascade states 
S^, particles is raised, = A + fig- All other excitation energies are independent 
of Us- As a consequence we find gapless {p,n) and (S~,S°)~^ excitations at fig = A. 
The situation is more complicated when kaon condensation is taken into account. In 
the kaon condensed phase there is mixing in the {p, S+, S~, S~) and {n, A*^, A°) 

sector. For m% -C /is -C A the spectrum is given by 



r, 







doN + i[To,N], 

{e (do + iXr) + (9o + iXl) e} 



(18) 




(19) 
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Figure 3: Left panel: Energy density as a function of the current Jk for several 
different values of fig = t^1I{'^Pf) close to the phase transition. Right panel: Ground 
state energy density as a function of /i^. We show the CFL phase, the kaon condensed 
CFL (CFLK) phase, and the supercurrent state (curCFLK). 

Numerical results for the eigenvalues are shown in Fig. El We observe that mixing 
within the charged and neutral baryon sectors leads to level repulsion. There are 
two modes that become light in the CFL window /i^ < 2A. One mode is a linear 
combination of proton and S"^ particles, as well as H~ and holes, and the other 
mode is a linear combination of the neutral baryons (n, 'EP , A^, A°). 

3.3 Meson supercurrent state 

What happens when gapless fermions appear in the spectrum? Several authors have 
shown that gapless fermion modes lead to instabilities in the current-current correla- 
tion function pO| 1^ . Motivated by these results we have examined the stability of 
the kaon condensed phase against the formation of a non-zero current |31 S] • Consider 
a spatially varying f/(l)y rotation of the kaon condensate 

/I \ 

U{x)iK^^U\x) = 1/V2 ze^<^^(^)/v^ . (20) 
V ie-^<^^(^VV2 1/V2 J 

This state is characterized by non-zero currents V and A. In order to determine 
the stability of the CFLK state we compute the vacuum energy as a function of the 
kaon current jk = V(j)K- The meson part of the effective lagrangian gives a positive 
contribution 

£ = \vlfjK. (21) 

A negative contribution can arise from gapless fermions. In order to determine this 
contribution we have to calculate the fermion spectrum in the presence of a non-zero 
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current. The spectrum is determined by the effective lagrangian f|T6l) . The dispersion 
relation of the lowest mode is approximately given by 

u;i = A + - -/i, - -V ■ 3K. (22) 

where I = v ■ p — pp and we have expanded ui near its minimum Iq = [fis + v ■ 
Jk)/4:. Equation (1221) shows that there is a gapless mode if fig > 4A/3 — jx/3. The 
contribution of the gapless mode to the vacuum energy is 

£ = ^JdlJ^u,9{-u;,), (23) 

where dQ is an integral over the Fermi surface. In Fig. [3] we show the ground state 
energy as a function of the current. We observe that there is a first order transition 
to a state with a non-zero current. The ground state energy as a function of Hs is 
shown in Fig. |3l see |l2l |13] for more details. Once the current becomes large the 
effective theory ceases to be reliable and states with multiple currents may appear. 
These states can be thought of as continuously connected to the crystalline quark 
matter phase fB]. 

4 Transport Properties 

Non-equilibrium properties, such as shear viscosity, bulk viscosity, thermal conduc- 
tivity and neutrino emissivity, play an important role in constraining the structure 
of compact stars. Shear and bulk viscosity control r-mode instabilities which, if not 
suppressed by viscous damping, lead to a fast spin-down of rapidly rotating compact 
stars. Neutrino emissivity controls the cooling behavior of the star, and neutrino 
opacities determine the spectral shape of the initial neutrino burst. In addition to 
these specific constraints there has been significant progress, in both theory and ob- 
servation, of tying together the rotation of the star, the magnetic field, and thermal 
properties. Exploiting these connections in order to constrain the phase structure of 
compact star matter will require a detailed understanding of transport properties. 

4.1 Hydrodynamics of the CFL phase 

The spontaneous breaking of ^/(l)^ is related to superfluidity, and the U{1)b effective 
theory can be interpreted as superfluid hydrodynamics [32]. We can define the fluid 
velocity as 

Va = D^ip, (24) 

where Da^f = + (/i, 0) and /io = {Da^D°'Lpy^'^. Note that this definition ensures 
that the flow is irrotational, V x t7 = 0. The identification flMl) is motivated by the 
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fact that the equation of motion for the U{1) field if can be written as a continuity 
equation 

d'^inova) = 0, (25) 

where no = S/iQ/vr^ is the superfiuid number density. At T = the superfiuid density 
is equal to the total density of the system, n = dP/dfi\^=f^Q. The energy-momentum 
tensor has the ideal fluid form 

Tap = (e + P)VaVl3 " PgafS, (26) 

and the conservation law d'^T^p = corresponds to the relativistic Euler equation 
of ideal fluid dynamics. We conclude that the effective theory for the U{1)b Gold- 
stone mode accounts for the defining characteristics of a superfiuid: irrotational, 
non-dissipative hydrodynamic flow. 

At non-zero temperature the hydrodynamic description of a superfiuid contains 
dissipative terms. Similar to the two fluid model of liquid helium we can describe the 
CFL phase, or any other relativistic superfiuid, as a mixture of an ideal superfiuid 
and a dissipative normal component [HI HS]. We will denote the densities of the 
superfiuid and normal components by ps and p„. We also define to be the velocity 
of the normal component, and to be the difference between the superfiuid and 
normal velocities. The normal fluid provides both non-dissipative and dissipative 
contributions to the energy momentum tensor. In the rest frame of the normal fluid 
the dissipative terms are 

ST, = {d.u, + d,u. - Is.Au,) - 5,CA iPsW,) - 5,C2d,u, (27) 

and 

6Toi = -KdiT . (28) 

Here t] is the shear viscosity, C,i^2 are bulk viscosities, and n is the thermal conductivity. 
Two additional bulk viscosities, ^3,45 control dissipative corrections to the equation of 
motion for the superfiuid velocity. There is a symmetry relation between the kinetic 
coefficients that requires that ^4 = Ci- Note that in the normal phase there is only 
one bulk viscosity, C, = (,2- 

4.2 Transport Coefficients 

The normal fluid is composed of quasi-particle excitations. In the CFL phase all quark 
modes are gapped and the relevant excitations are Goldstone bosons. At very low 
temperature, transport properties are dominated by the massless Goldstone boson 
associated with the breaking of the U{1)b symmetry. The effective lagrangian ([6]) 
determines the rates for the relevant scattering processes. 
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Figure 4: Bulk viscosity ^ = ^2 as a function of temperature for an oscillation period 
T = 2'k/u = 1 ms. CFL phase: contribution from the process v? + 93 for 

different values of 5m = 171^0 — fis and contribution from ip -(r^ ip + ip. 2SC phase and 
unpaired quark matter: contribution from the process u + d -H- u + s. 

Shear viscosity is related to momentum transport. At low temperature the shear 
viscosity of the CFL phase is determined hjip + ip^ip + ip scattering. Manuel et 
al. find US] 

^^1.3x10-^^. (29) 

The bulk viscosity is sensitive to particle number changing processes. This includes 
purely strong decays like (p ^ ip + ip, or electroweak processes like the strangeness 
changing reaction Lp + ip. We first consider the pure QCD contribution. Bulk 

viscosity vanishes in an exactly scale invariant system. For realistic quark masses the 
dominant source of scale breaking is the strange quark mass. The contribution from 
the process (p -H- ip + ip is [17] 

4 

(2^0.011^. (30) 

We show this contribution in Fig. HI The electroweak process Lp + Lp was 

studied in [18]. The weak contribution has a significant frequency dependence. In 
Fig. Hjwe show the results for an oscillation period t = 27t/u = 1 ms. We observe 
that at T ~ (1 — 10) MeV the bulk viscosity of CFL matter is comparable to that of 
unpaired quark matter. For T < 1 MeV, (2 is strongly suppressed. Depending on the 
poorly known value for Sm = ttiko — /is the pure ip contribution given in equ. f l30|) 
may dominate over the ^ if + ip reaction at low enough temperatures. However, 
for T < 0.1 MeV the ip mean free path is on the order of the size of the star, i.e., the 
system is in the coUisionless rather than in the hydrodynamic regime, and the result 
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ceases to be meaningful. 

The thermal conductivity of a CFL superfluid was studied by Braby et al. |19] . The 
calculation is subtle because k vanishes for a system of quasi-particles with exactly 
linear dispersion relations |50j. The reason is that k measures the rate of energy 
transport relative to the motion of the fluid, but in a gas of massless particles with 
linear dispersion one cannot transport energy without transporting momentum. As 
a consequence, thermal conductivity is sensitive to non-linearities in the dispersion 
relation. Braby et al. find |l9] 

K ~ 4.01 X 10^2 Me (31) 

They also estimate the contribution to /t from phonon scattering on kaons. This term 
grows as a/T, but it is significantly smaller than the phonon contribution in the regime 
where the calculation is reliable. Non-linearities in the dispersion relation also play 
a role in determining the remaining two bulk viscosities, C,i and Mannarelli 
and Manuel find Ci ~ f^lKl^T) and ^3 ~ l/(T/i^). Note that is non-zero even in 
the approximately conformal limit — )■ 0. 

4.3 Neutrino emissivity 

In CFL quark matter the neutrino emissivity is dominated by reactions involving 
pseudo-Goldstone modes such as 



e ) 



7r° ^ Z/e+I7e, (32) 
(p + (f (p + Ue+Ve. 

These processes were studied in [52| [53] . The decay rates of the massive mesons vr^, 
K"^, and 7r° are proportional to their number densities and are suppressed by Boltz- 
mann factors exp{—E/T), where E is the energy gap of the meson. The emissivity 
from ir"^ decay is 

' iGlf.ml)mlnJl + 2(l-vl) + ^^(l-vlY], (33) 



in ' ' ' \ ^ vim? 



where n,r is the number density of pions. Similar results can be derived for 7r° and 
decay. Since the pseudo-Goldstone boson energy gaps are on the order of a few 
MeV, the emissivities are strongly suppressed as compared to unpaired quark matter 
for temperatures below this scale. Neutrino emission from processes involving the Lp 
is not exponentially suppressed, but it involves a very large power of T, 

--%r. (34) 
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and is numerically very small. Reddy et al. also studied the neutrino mean free path 
lu. For T ~ 30 MeV the mean free path is on the order of 1 m, but for T < 1 MeV, 
> 10 km ^ 



5 Outlook 

There are a variety of issues that remain to be studied. While the calculation of 
transport coefficients in the CFL phase is now essentially complete, this is not the 
case for many of the less dense phases. There are calculations of shear and bulk 
viscosity as well as neutrino emissivity in the CFL-K phase [HU [551 ES] , but there are 
essentially no results for the spatially inhomogeneous or anisotropic phases. There is 
also much work to be done in order to understand many phenomena that are relevant 
in compact stars, like mutual friction between the normal fluid and superfluid vortices 
\57\ EH] , or the role of the quark-hadron interface. 
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